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ABSTRACT 

We report the first detection of pulsed X-ray emission from the young, energetic radio and 7-ray 
pulsar PSR B1706-44. We find a periodic signal at a frequency of / = 9.7588088 ± 0.0000026 Hz 
(at epoch 51585.34104 MJD), consistent with the radio ephemeris, using data obtained withthe High 
Resolution Camera on-board the Chandra X-ray Observatory. The probability that this detection is a 
chance occurrence is 3.5 x 10 -5 as judged by the Rayleigh test. The folded light curve has a broad, 
single-peaked profile with a pulsed fraction of 23% ± 6%. This result is consistent the ROSAT PSPC 
upper limit of < 18% after allowing for the ability of Chandra to resolve the pulsar from a surrounding 
synchrotron nebula. We also fitted Chandra spectroscopic data on PSR B1706— 44, which require at 
least two components, e.g., a blackbody of = (1.66^q Jg) x 10 6 K and a power-law of V = 2.0 ± 0.5. 
The blackbody radius at the nominal 2.5 kpc distance is only = 3.6 ± 0.9 km, indicating either a 
hot region on a cooler surface, or the need for a realistic atmosphere model that would allow a lower 
temperature and larger area. Because the power-law and blackbody spectra each contribute more than 
23% of the observed flux, it is not possible to decide which component is responsible for the modulation 
in the spectrally unresolved light curve. 

Subject headings: pulsars: individual (PSR B1706— 44) — stars: neutron — X-rays: stars 



age and spin-down luminosity. Finley et al. (1998), using 
data from the ASCA CIS, placed an upper limit of 22% on 
the pulsed fraction in the 2—10 keV band, although their 
image was severely contaminated by stray light from the 
bright low-mass X-ray binary 4U1705— 44 just outside the 
field of view. Finley et al. also analyzed a ROSAT HRI 
observation, from which they concluded that 57 ± 12% of 
the photons associated with PSR B1706— 44 actually come 
from a compact nebula with an exponential scale length of 
~ 27". When folding only the photons from the HRI point 
source at the radio period, they derived an upper limit of 
29% on its pulsed fraction. An observation by the non- 
imaging Rossi X-ray Timing Explorer performed during a 
low state of 4U1705— 44 also failed to detect pulsations in 
the 9-18.5 keV band. (Ray, Harding, & Strickman 1999). 
Here we report on Chandra observations of PSR B1706— 44 
in which we detect its pulsed emission in X-rays. 

2. OBSERVATIONS 

The field of PSR B1706-44/G343.1-2.3 was observed 
by Chandra twice, once with each of the two types of cam- 
eras at the focal plane of the telescope. Herein we analyze 
data from these observations, made available though the 
Chandra public archive. 

Timing data on PSR B1706— 44 were acquired on 2000 
Feb 11 with the imaging High Resolution Camera (HRC-I; 
Murray et al. 1997), a multichannel plate detector with 
an effective time resolution degraded to ~ 4 ms due to a 
known timing error. The HRC is sensitive to X-rays in 
the energy range 0.1 — 10 keV; no useful spectral infor- 
mation is available. The target was placed near (« 28") 
the optical axis where the telescope mirror point-spread 
function (PSF) has half power radius (the radius enclos- 
ing 50% of the total source counts) of ~ 0."5 for energies 
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1. INTRODUCTION 

PSR B1706— 44 is a 102 ms radio pulsar discov- 
ered by Johnston et al. (1992) and possibly associ- 
ated with the supernova remnant G343. 1—2.3 (McAdam, 
Osborne, & Parkinson 1993). With characteristic age 
P/2P — 17,500 yr and spin-down luminosity E = 3.4 x 
10 36 ergs s _1 PSR B1706— 44 is considered a young, ener- 
getic neutron star. The distance to the pulsar is some- 
whate uncertain. According to the Taylor & Cordes 
(1993) free electron model of the Galaxy the pulsar lies 
1.8 kpc away, however Koribalski et al. (1995) find a 
kinematic distance in the range 2.4 — 3.2 kpc from H I 
absorption. Originally detected as the high-energy 7-ray 
source 2CG342-02 by the COS B satellite (Swanenburg 
et al. 1981), PSR B1706-44 is one of approximately eight 
rotation-powered pulsars that are responsible for some of 
the brightest 7-ray sources in the sky. Thompson et al. 
(1992; 1996) showed that photons detected from the di- 
rection of PSR B1706-44 by the EGRET instrument on 
board the Compton Gamma-Ray Observatory are pulsed 
at the radio period. All of the known 7-ray pulsars are well 
observed at X-ray energies between 0.1 and 10 keV, where 
both surface thermal emission and magnetospheric syn- 
chrotron emission can be studied. However, pulsed X-ray 
emission has not been detected from PSR B1706— 44 de- 
spite several searches. This failure has hindered attempts 
to understand its X-ray emission mechanism(s). 

Using an observation of PSR B1706-44 by the ROSAT 
PSPC, Becker, Brazier, Trumper (1995) placed an upper 
limit of 18% on the pulsed fraction. They attributed this 
negative result to the diluting effect of an as-yet unre- 
solved synchrotron nebula, recalling the history of Vela 
(Ogelman, Finley, & Zimmerman 1993), a pulsar of similar 
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E < 6 keV. The HRC-I oversamples the PSF with pixels 
0."1318 on a side, which allows the pulsar to be isolated 
from any surrounding emission such as the pulsar wind 
nebulae (PWNe) typically found associated with young 
pulsars observed by Chandra (e.g., see Gotthelf 2001). 

Data were also obtained with the Advanced CCD Imag- 
ing Spectrometer (ACIS; Burke et al. 1997), which has 
resolution AE/E ~ 0.1 at 1 keV scaling as l/VE over its 
0.2-10 keV active band-pass. The pulsar was observed on 
2000 Aug 14 and positioned on the back-illuminated S3 
chip of the ACIS-S array, offset by 0'5 from the aim-point, 
where the PSF is undersampled by the 0."4920 x 0."4920 
CCD pixels. Data were collected in the nominal timing 
mode, with 3.241 s exposures between CCD readouts, and 
in "FAINT" spectral mode. 

The standard Chandra screening criteria produced a to- 
tal usable exposure time of 46.8 ks and 14.3 ks for the 
HRC and ACIS data sets, respectively. These resulted in 
a similar number of detected photons from the source, and 
the images made from the two instruments are found to 
be consistent. The ACIS image, restricted to the 0.3 — 8 
keV energy band to reduce instrumental background, re- 
veals a point-like X-ray source at the radio pulsar posi- 
tion and four nearby faint X-ray sources that are coin- 
cident with USNO stars to sub-arcsecond accuracy. The 
Chandra position of PSR B1706-44 is a = 17 h 09 m 42 s .73, 
S = -44°29'08'.'4 (J2000.0) with RMS error 0"5. This is in 
agreement with the precise radio timing position of Wang 
et al. (2000, see Table 1). 

To look for evidence of extended emission around 
PSR B1706— 44 we constructed a radial profile, the dis- 
tribution of counts per unit area in the HRC as a function 
of distance from the pulsar. For comparison, we analyzed a 
23 ks Chandra observation of the millisecond pulsar PSR 
J0437-4715 observed at the aim-point of the HRC just 
two days later (Zavlin et al. 2002). PSR J0437-4715 is 
an isolated point source with a spectrum similar to that 
of PSR B1706— 44 (as shown in §4), thus providing a re- 
alistic example of the in-orbit point-spread function. Fig- 
ure |l] displays the two radial profiles after adjusting for 
different background intensities, sampled in an annulus of 
30" < r < 60" where the background is flat. We then 
normalized the two profiles to the peak intensity within 
r < l."0. While the majority of the emission is point- 
like, there is clearly diffuse emission out to a radius of 

20". There is also a bump of enhanced emission be- 
tween 11" < r < 23". This result is consistent with the 
analysis of these data by Dodson & Golap (2002), who dis- 
played extended emission out to w 5", and is marginally 
consistent with the decomposition of the ROSAT HRI im- 
age deduced by Finley et al. (1998). In order to make 
a direct comparison with the ROSAT energy band, we 
also extracted photons from the ACIS image in the range 
0.2 — 2.0 keV. After background subtraction, we find 513 
photons in the ACIS from the point source (radius l."5), 
and 383 nebular photons in the annulus l."5 < r < 20". 

3. TIMING ANALYSIS 

We searched for a pulsed signal from PSR B1706— 44 in 
the HRC data using time-tagged photon events extracted 
from an 3."0 radius aperture centered on the X-ray pulsar 
position. This aperture effectively excludes most of the dif- 



fuse emission from the putative surrounding synchrotron 
nebula; furthermore, less than 3% of the photons in this 
aperture are background events. The arrival times of the 
824 selected photons were corrected to the solar system 
baryccntcr using the JPL-DE200 cphcmcris and the radio 
timing position. We then generated a periodogram using 
the Z\ statistic (Rayleigh test) on a range of test frequen- 
cies centered on the predicted pulsar frequency derived 
from a fit to radio observations from the Parkes Obser- 
vatory archives (Wang et al. 2000) that span the X-ray 
observation. The residual deviations of the radio pulse ar- 
rival times from the fitted quadratic ephemeris in Table 1 
are less than 1 ms, including one observation that occurred 
24 days before the X-ray observation. The known disper- 
sion measure was used to correct the radio pulse arrival 
times to infinite frequency for absolute phase comparison 
with the X-ray pulse. 

The X-ray periodogram is shown in Figure |[ We 
find a peak at / = 9.7588088 ± 0.0000026 Hz (epoch 
51585.34104 MJD) which, within the quoted 68% uncer- 
tainty range, is identical to the contemporaneous radio fre- 
quency, / = 9.7588097 Hz. The Z\ statistic for this peak 
is 20.54, which has a probability of chance occurrence of 
3.5 x 10~ 5 . The light curve, shown in Figure 0, is broad 
and single peaked with a pulsed fraction of 23%± 6%. Ac- 
cording to convention, the pulsed fraction is defined as the 
ratio of counts above the minimum in the light curve to 
the total counts. This detection is consistent with the pre- 
vious upper limits of < 29% from the ROSAT HRI (Finley 
et al. 1998) and < 18% from the ROSAT PSPC (Becker 
et al. 1995). The latter did not resolve the pulsar from 
the nebula, therefore, it represents a true upper limit of 
< 31% given the contribution of diffuse photons measured 
in the Chandra ACIS image as described in §2. The arrival 
time of the radio pulse is consistent with occurring at the 
center of the broad X-ray peak. This is not the case for 
the EGRET 7-ray pulse, which is centered « 0.37 cycles 
after the radio pulse (Thompson et al. 1996). 

4. SPECTROSCOPY 

To characterize the energy dependence of the emission 
from the pulsar and its putative wind nebula, we analyzed 
spectral data obtained with the ACIS detector. Pulsar 
and PWN source counts spectra were extracted from two 
concentric regions, a circle of radius 1" and an annulus 
of radii 2" < r < 10", respectively. An annular back- 
ground region centered on the pulsar was extracted from 
radii 30" < r < 60", which appears to be pure background 
(see Fig. [l]). We verified that background is a negligible 
contribution to the pulsar spectrum and a 10% contribu- 
tion to the PWN spectrum over their respective extrac- 
tion regions. To be specific, there are 561(486) source 
counts and an estimated 1(54) background count (s) for 
the pulsar(PWN) regions in the 0.5 — 8.0 keV spectral fit- 
ting range. From the radial profile (Fig. |l|) it is clear that 
there is negligible cross-contamination between the pulsar 
and PWN extraction regions. 

The ACIS pulse-height data were corrected for gain and 
resolution degradation due to CCD charge transfer ineffi- 
ciencies (CTI) using the correctit software (Townsley et al. 
2001). We used the supplied instrument response matrices 
appropriate for the CTI corrected data. Since the loca- 
tion of the target on the CCD straddled regions covered 
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by two response files, we generated a count-weighted mean 
response to match the data. The CIAO tool mkarf were 
used to generate a point-source mirror response matrix for 
the pulsar and PWN source regions. The extracted pho- 
tons were binned in pulse-height channel space to match 
the response matrices and regrouped such that each fit- 
ted spectral bin contained a minimum of 20 counts. The 
resulting spectra were then fitted with the X-ray spectral 
analysis package XSPEC (version 11.1) using three different 
models, a power law, a simple blackbody, and a sum of 
both, each with interstellar absorption. 

Initial fits to the pulsar spectrum using blackbody or 
power-law models were unacceptable, with % 2 = 115 and 
X 2 — 60 for 23 degrees-of-freedom (DoF), respectively. We 
next determined, as is often the case for young pulsars, 
that a blackbody plus power-law fit provides an adequate 
description of the spectrum. The intrinsic parameters of 
this multicomponent model are rather unconstrained and 
strongly correlated with the fitted value of the interstellar 
absorption column density, Nh ■ To remove an unnecessary 
degree of freedom, we fixed iVn at the value obtained from 
a simple fit to the PWN, 5.5 x 10 21 cm" 2 , since we expect 
the intervening column density to be the same. The PWN 
is well characterized by a single absorbed power-law with 
X 2 = 13 for 20 DoF. The results of these fits are given in 
Table and the counts spectrum of the pulsar, with the 
model components superposed, are shown in Figure ^. 

The fitted blackbody component yields effective tem- 
perature Too = (1.66jZ ^ 15 ) x 10 6 K, bolometric luminos- 
ity Loo = 6.81°;® x 10 32 , and effective radius = 
3.6 ± 0.9 km in the observed frame. Here we adopt a 
distance of 2500 pc as a compromise between the smaller 
DM distance of 1.8 kpc, and the H 1 kinematic distance of 
2.4—3.2 kpc. While Too is consistent with some standard 
neutron star cooling curves (Umcda, Tsuruta, & Nomoto 
1994), the associated i?oo is significantly smaller than the- 
oretical neutron star radii. As has been found frequently 
for other pulsars, e.g., Vela (Pavlov et al. 2001), such an 
outcome could have either of two implications. First, the 
thermal X-rays could be coming from a smaller hot region 
on the surface of a cooler neutron star, the substantial in- 
terstellar iVn making it difficult to measure softer X-rays 
coming from the full surface. Second, a more realistic at- 
mosphere model, especially if dominated by hydrogen or 
helium, would require a lower effective temperature and 
larger effective radius, perhaps becoming consistent with 
emission from the full neutron star surface. 

5. DISCUSSION AND CONCLUSIONS 

There are at least two mechanisms for the emission of 
broad-band pulsed X-rays from young rotation-powered 
pulsars like PSR B1706— 44. One is non-thermal mag- 
netospheric synchrotron from relativistic electrons and 
positrons created either in regions above the neutron star 
polar caps or in outer gaps. The second is thermal emis- 
sion from the hot surface, a result of initial cooling of the 
hot neutron star or reheating of the polar caps by back- 
flowing accelerated particles. Sometimes the shape of the 



pulse sheds additional light on the X-ray emission mech- 
anism. Sharp, narrow pulses of high amplitude can only 
be produced by a highly beamed, thus relativistic popula- 
tion of electrons, while quasi-sinusoidal pulses of low am- 
plitude such as describe PSR B1706— 44 can be produced 
by either mechanism. For those intermediate-age pulsars 
(10 4 — 10 6 yr) that are also EGRET sources, the pres- 
ence of both types of X-ray source is usually discovered 
when spectrally resolved timing data are available (e.g., 
Wang et al. 1998; Pavlov et al. 2001). Unfortunately, in 
this case, the Chandra HRC has little or no energy reso- 
lution. Since each spectral component fitted to the ACIS 
spectrum contributes more than 23% of the flux in the to- 
tal energy band to which the HRC is sensitive (60% from 
blackbody, 40% from power law), the pulsed fraction alone 
does not reveal the source of the pulsed X-rays. Either 
or both components may contribute to the modulation. 
Future spectrally resolved X-ray observations with high 
throughput and moderate time resolution, such as with 
XMM-Newton, could resolving this ambiguity. 

The > 100 MeV luminosity of most 7-ray pulsars is 
a significant fraction of their spin-down power. In the 
case of PSR B1706-44 this fraction is w 0.20 if isotropic, 
while its 0.5 — 8 keV non-thermal X-ray luminosity is only 
1.3 x 10~ 4 E including both pulsar and nebula, similar to 
that of other pulsars. Thompson et al. (1996) parame- 
terized the EGRET spectrum of PSR B1706-44 as a bro- 
ken power law, with photon index Y = 1.27 ± 0.09 from 
50 MeV to 1 GeV, steepening to V = 2.25 ± 0.13 above 
1 GeV. Since there is no evidence for any unpulsed 7-ray 
emission, the EGRET spectrum can be compared directly 
with the power-law component of the pulsar point source 
in the Chandra ACIS spectrum. When extrapolated back 
to 10 keV, the EGRET spectrum nearly matches the X- 
ray flux, although the power-law X-ray slope itself, with 
r = 2.0±0.5, is somewhat steeper than the EGRET value. 
The non-thermal X-rays are not likely to be emitted by 
the same population of electrons/positrons as produce the 
7-rays, but may instead originate from a much less ener- 
getic population, such as secondary pairs from the inward 
emitted 7-rays converting on the strong B-field near the 
neutron star surface. This distinction is also supported 
by the observed phase offset between the X-ray and 7-ray 
pulses (§3). The Wang et al. (1998) model predicts L x 
(2-10 keV) w 2 x 10 31 ergs s _1 and T = 1.5 from this pro- 
cess for PSR B1706— 44, which is similar to its observed 
non-thermal X-ray component. If its thermal luminosity 
comes partly from a small surface area, then inward flow- 
ing primary electrons impacting the polar caps may be 
responsible. Wang et al. predicted L x ^ 1 x 10 33 ergs s~ 4 
of thermal emission from this process, not far from the 
observered thermal luminosity of PSR B1706— 44. 



This work made use of data obtained from the Chandra 
Data Archive, and is supported by NASA LTSA grant 
NAG 5-7935 to EVG. We thank D. Lewis and R. N. 
Manchester for their assistance with the radio ephemeris. 



4 



Pulsed X-ray Emission from PSR B1706-44 



REFERENCES 



Becker, W., Brazier, K. T. S., & Triimper, J. 1997, A&A, 298, 528 
Burke, B. E., Gregory, J., Bautz, M. W., Prigozhin, G. Y., Kissel, 

S. E., Kosicki, B. N., Loomis, A. H., & Young, D. J. 1997, IEEE 

Transac. Elec. Devices, 44, 1633 
Dodson, R., & Golap, K. 2002, MNRAS, submitted 
Finley, J. P., Srinivasan, R. , Saito, Y., Hirayama, M., Kamae, T., &: 

Yoshida, K. 1998, Ap. J., 493, 884 
Gotthelf, E. V. 2001, in "Proceedings pf the Texas Rfjlativistic 

Astrophysics Meeting," Austin, Texas ( |astro-ph/010512i ) 
Johnston, S., Lyne, A. G., Manchester, R. N., Kniffen, D. A., 

D'Amico, N., Lim, J., & Ashworth, M. 1992, MNRAS, 255, 401 
Kaspi, V. M., Bailes, M., Manchester, R. N., Stappers, B. W., 

Sandhu, J. S., Navarro, J., & D'Amico, N. 1997, ApJ, 485, 820 
Koribalski, B., Johnston, S., Weisberg, J. H., & Wilson, W. 1995, 

ApJ, 441, 756 

McAdam, W. B., Osborne, J. L., & Parkinson, M. L. 1993, Nature, 
361, 516 

Murray, S. S., et al., 1997, SPIE, 3114, 11 



Ogelman, H., Finley, J. P., & Zimmerman, H. U. 1993, Nature, 361, 
136 

Pavlov. G. P., Zavlin, V. E., Sanwal, D., Burwitz, V., & Garmire, 

G. P. 2001, ApJ, 552, L129 
Ray, A., Harding, A. K., & Strikman, M. 1999, ApJ, 513, 919 
Swanenburg, B. N., et al. 1981, ApJ, 243, L69 
Taylor, J. H., & Cordes, J. M. 1993, ApJ, 411, 674 
Thompson, D. J., et al. 1992, Nature, 359, 615 
Thompson, D. J., et al. 1996, ApJ, 465, 385 

Townsley, L. K., Broos, P. S., Garmire, G. P. & Nousek, J. A. 2000, 

ApJ, 534, L139 
Umeda, H., Tsuruta, S., & Nomoto, K. 1994, ApJ, 433, 256 
Wang, F. Y.-H., Ruderman, M., Halpern, J. P., & Zhu, T. 1998, 

ApJ, 498, 373 

Wang, N, Manchester, R. N., Pace, R. T., Bailes, M., Kaspi, V. M., 

Stappers, B. W., & Lyne, A. G. 2000, MNRAS, 317, 843 
Zavlin, V. E., Pavlov, G. G., Sanwal. D., Manchester, R. N, 
Mmpor, J M Halp ern, J. P., Sz Becker, W. 2002, ApJ, 000, L000 



( astro- ph/01 12544) 



Gotthelf et ai 



5 



Table 1 

Radio Ephemeris of PSR B1706-44 



Parameter 



Value 



R.A. (J2000) 

Decl. (J2000) 

Valid Range (MJD) 

Epoch (MJD) a 

Frequency (Hz) 

Frequency Derivative (Hz s -1 ) 

Dispersion Measure (cm~ 3 pc) 



17 h 09 m 42?728 
-44°29'08."24 
51520- 51956 
51585.3410392671 
9.7588096779 
-8.880796 x 10~ 12 
75.69 



a TDB epoch of radio peak at infinite frequency. 



Table 2 

Spectral Fits to PSR B1706-44 and its PWN a 



Pulsar Nebula 



JV H (10 21 cm" 2 ) 


5.5 (fixed) 


5 5+ 2 - 5 

°-°-2.1 


r 


2 +0 - 5 

z - u -0.5 


1 Q4+0.24 
i '°*-0.30 


C p i h 


4.3±f- g x 10~ 5 


5.5±?-g x 10" 5 


F p i (absorbed) 


1.2 x 'lO" 13 


3.4 x 10" 13 


L p i (unabsorbed) d 


1.45±g;$ x 10 32 


3.lt\i x 10 32 


Toe (K) c 


lMtall x 106 




Roc (km) c 


3.6 ± 0.9 




Too (ergs s" 1 ) 


6.8±?"? x 10 32 




Xl [DoF] 


1.4 [21] 


0.61 [21] 



a Uncertainties are 68% confidence for two interesting parameters. 
b Powcr-law normalization in units of photons cm~ 2 s _1 kcV -1 at 1 keV. 
c Power-law flux (ergs cm -2 s _1 ) in the 0.5-8 keV band. 
d Powcr-law luminosity (ergs s — 1 ) in the 0.5-8 keV band for d = 2500 pc. 
°Sphcrical blackbody parameters for d = 2500 pc. 
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Fig. 1. — X-ray intensity of PSR B1706— 44 as a function of radius 
measured with the Chandra HRC (histogram), compared with the 
HRC response to a point source (solid line; see text). The radial 
profile shows clear evidence for extended emission between l."5 and 
20" . The error bars are la. 
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Fig. 3.— The light curve of PSR B1706-44 from the Chandra 
HRC folded according to the radio ephemeris of Table 1. The quasi- 
sinusoidal modulation has pulsed fraction of 23% ± 6%. Phase zero 
corresponds to the epoch of the radio peak. The error bars are la. 
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Fig. 2. — Evidence for pulsed X-ray emission from PSR B1706— 44 
using a periodogram. The arrow indicates the frequency expected 
according to the extrapolated radio ephemeris (see text). The inset 
is an expanded view around the peak. The peak value = 20.54 
at / = 9.7588088 Hz has a 3.5 X 10~ 5 probability of occurring by 
chance. 
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Fig. 4. — Chandra ACIS-S3 spectrum of the point source coincident 
with PSR B1706— 44. Top panel: Data (crosses) and best-fit model 
(dark line) for the parameters given in Table be thin lines show the 
contribution of the power-law (PL) and blackbody (BB) components 
of the fit. Bottom panel: Difference between the data and model, in 
the same units as the top panel. 



